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ABSTRACT 

Stellar evolution models predict the existence of hybrid white dwarfs (WDs) with a carbon- 
oxygen core surrounded by an oxygen-neon mantle. Being born with masses ^1.1 Mq, hy¬ 
brid WDs in a binary system may easily approach the Chandrasekhar mass (Mch) by accretion 
and give rise to a thermonuclear explosion. Here, we investigate an off-centre deflagration in 
a near-Mch hybrid WD under the assumption that nuclear burning only occurs in carbon-rich 
material. Performing hydrodynamics simulations of the explosion and detailed nucleosynthe¬ 
sis post-processing calculations, we find that only 0.014 Mq of material is ejected while the 
remainder of the mass stays bound. The ejecta consist predominantly of iron-group elements, 
O, C, Si and S. We also calculate synthetic observables for our model and find reasonable 
agreement with the faint Type lax SN 2008ha. This shows for the first time that deflagrations 
in near-Mch WDs can in principle explain the observed diversity of Type lax supernovae. 
Leaving behind a near-Mch bound remnant opens the possibility for recurrent explosions or 
a subsequent accretion-induced collapse in faint Type lax SNe, if further accretion episodes 
occur. From binary population synthesis calculations, we And the rate of hybrid WDs ap¬ 
proaching Mch to be on the order of 1 percent of the Galactic SN la rate. 

Key words: supernovae: individual: SN 2008ha - methods: numerical - hydrodynamics - 
radiative transfer - nuclear reactions, nucleosynthesis, abundances - stars:evolution 


1 INTRODUCTION 

Type lax supemovae (SNe lax) form a distinct class of astronom- 
ical transients ([Foley et al.|2013|> w ith SN 2002cx jLi et al.|2003[ 
[Branch et al.|2004 Jha et al.|2006| l as the proto-typical example. It 
has been argued that SNe lax occur at a rate of 5 to 30 percent of the 
overall rate of Type la supemovae (SNe la) (|Li et al.|201 lb||Foley| 
|et al.|2013[ [White et al.|20T5| >. Recently, pdcCully et al | ( |2014^ re- 
ported the first detection of a progenitor candidate for a SN lax. In 
pre-explosion Hubble Space Telescope (HST) images they discov¬ 
ered a point source at the location of SN 2012Z. Future HST ob¬ 
servations will show whether this point source is actually a massive 
star that exploded giving rise to the SN, or a He donor to an accret- 
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ing carbon-oxygen (CO) white dwarf (WD) as argued by|McCully| 
[etaL] ( |20T4) . 

Although spectroscopically similar to SNe la, SNe lax are 
characterized by distinctive features: (i) expansion velocities at 
maximum light (2000 to 8000 km s“^) that are significantly lower 


than those of SNe la (10,000 to 15,000kms“^, e.g. [Foley et al. 
[2013[ [Stritzinger et al.[[2015[ l, (ii) low peak absolute magnitudes 
with respect to the width-luminosity relation of SN la light curves 
(e.g. [ McClelland et al.|2Q10[[Stritzinger et al.|2014[ l, (iii) peculiar 
late-time spectra that differ strongly from other SNe (e.g. [Jha et al. [ 
[2006[[S^u et al.|2008| > and (iv) a statistical association with star- 
formation regions (so far no SN lax has been observed in elliptical 
galaxies), indicating short delay times of 30 to 50 Myr ([Lyman et al.j 
|20T3l (. 

SNe lax show significant diversity in peak absolute magni- 
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tude and expansion velocities. For example, SN 2005hk, one of the 
brighter class members, peaked at ~ —18.1 mag i Chornock 


|et al.|[200^ Phillips et al.|[2007{ |Sahu et ^|2008| l. In contrast, 
SN 2008ha, to date th e faintest member of the class, reached only 


m:- 


— 14.2 mag jpoley et ak 


2009 


Valenti et al. 


2009 


Foley 


|et al.|2010l l. This indicates a wide range of ejecta properties. 

Several explosion mechanisms have been proposed for 
SNe lax, including core collapse of massive stars jFoley et al.|2009| 
[Valenti et al.|2009[ |Moriya e t al.|2010| >, detonations in He shells on 
top of CO WDs I Foley et al.|2009| ) and deflagrations or pulsational 
delayed detonations in near Chandrasekhar-mass (Mch) CO WDs 
( [Branch et al.||2004^ [Phillips et al.[[200'7| [Stritzinger et al.[|2015^ . 

The mechanism that seems to agree best with the observed proper¬ 
ties of SNe lax, including the possible detection of a stellar remnant 
for the faint Type lax SN 2008ha ( [Foley et al.|2014b[ l, is a deflagra¬ 
tion in a near-Mch CO WD that fails to completely unbind the WD, 
leaving behind a bound remnant ( [Foley et al.|2009[ l. 

[Jordan et al.[ ( [2012t were the first to present hydrodynamical 
simulations of this scenario yielding ®®Ni and total ejecta masses 
in agreement with the observationally derived quantities for bright 
SNe lax. Combining detailed explosion models and radiative trans¬ 
fer simulations, [Kromer et al.| ( [2013^ have shown that pure defla¬ 
grations in near-Mch CO WDs that leave behind a compact rem¬ 
nant can indeed reproduce the observational display of the brighter 
SNe lax like SNe 2002cx or 2005hk. |Fink et al.|j2014[ l have con¬ 
firmed these results with a set of models exploring deflagrations in 
near-Mch CO WDs for different ignition setups. However, they fail 
to explain the faint class members like SN 2008ha (their faintest 
model yields = —16.8 mag while SN 2008ha peaks at 

— 14.2 mag; [Foley et al.|200^ . Although [Jordan et al.[p012} and 
[Fink et aT] 2014 1 did not fully explore the parameter space of 
possible ignition configurations, it seems implausible to reach low 
enough ®®Ni masses to explain the low luminosities of the faintest 
SNe lax from deflagrations in near-Mch CO WDs. 

To explain these events with deflagrations in near-Mch WDs, 
one somehow has to quench burning at an early stage to reduce the 
®®Ni production. This is potentially possible by an abrupt change in 
the composition profile of the exploding WD as seen in hybrid WDs 
where a CO core is surrounded by an oxygen-neon (ONe) man¬ 
tle. Such objects have recently been predicted by stellar evolution 
models jPenissenkov et al.|2013[|Chen et al.|2014[ l. |Denissenkov[ 
[et al.| ( (20T5l l have shown that a thermonuclear runaway in accreting 
hybrid CONe WDs is possible when they approach Mch- Since hy¬ 
brid CONe WDs originate from relatively massive zero-age main 
sequence (ZAMS) stars jChen et al.|2014[ l, such a scenario could 
also explain the short delay times of SNe lax ( [Men g & Podsiad^ 
[lowski|2014[ l. 

Here, we study a deflagration in such a near-Mch hybrid 
CONe WD. The paper is organised as follows: in Section]^ we 
discuss the chosen progenitor model before we present our explo¬ 
sion simulation in Section Synthetic observables for our model 
and comparison to observational data are presented in Section]^ In 
Sectionj^we discuss our results and constraints on rates and delay 
times from binary population synthesis calculations before present¬ 
ing our conclusion in Section]^ 


2 PROGENITOR MODEL 

In the single-degenerate progenitor scenario for SNe la, a CO WD 
is supposed to accrete H-rich material from a binary companion 
until it approaches Mch and a thermonuclear runaway occurs (e.g. 


[Hillebrandt & Niemeyer|2000[ for a review). CO WDs are formed 
when the cores of asymptotic giant branch (AGB) stars are exposed 
after they have expelled their envelopes. In classical single-star evo¬ 
lution models the maximum mass for the CO core of an AGB star 
is ~ 1.05 Mq ( [Chen et al.|2014[ l. For more massive cores, C burn¬ 
ing is ignited (super-AGB phase) and the entire core is converted 
to O and Ne, leading to the formation of an ONe WD once the star 
has lost its envelope. ONe WDs have been discarded as progenitor 
candidates for the single-degenerate progenitor scenario of SNe la. 
It has been argued that they should undergo a gravitational collapse 
and form a neutron star rather than a thermonuclear explosion when 


approaching Mch iNomoto|1984| 

1987| Nomoto & Kondo 

19911. 

Recently, Denissenkov et al. 

(2013J and Chen et al. 

]20T^ 


have shown that convective boundary mixing during the super- 
AGB phase can prevent off-centre ignited C burning from reaching 
the centre of the star, thus leaving a hybrid CONe WD with a C-rich 
core at the end of the super-AGB phase. In contrast to an ONe WD, 
such a hybrid WD with a C-rich core is likely to be able to trig¬ 
ger a thermonuclear explosion when approaching Mch- ONe-rich 
WDs containing a small fraction of unbumt C have been predicted 
by stellar evolution models previously ( [Siess|2006l l, and have been 
mentioned in the literature as possible progenitors of thermonuclear 
expl osions ([Garcia -Berro e t al.|1997t . 

[Denissenkov et al.[P015[ l have simulated the mass accretion 
onto hybrid CONe WDs to investigate their potential as SN la pro¬ 
genitors in the single-degenerate scenario. They find that a ther¬ 
monuclear runaway in accreting hybrid CONe WDs is possible de¬ 
pending on the mass of the parent star and accretion history. The ex¬ 
act thermal and chemical structure of the hybrid WD at the time of 
explosion, however, is uncertain and can only be determined from 
3D reactive-convective hydro simulations. In particular, a series of 
convective Urea shell flashes can either limit the mass of the con¬ 
vective C core or completely suppress convective mixing ( [Denis-[ 
[senkov et al.|2015) . 

Here, we investigate how a deflagration will incinerate a hy¬ 
brid CONe WD that undergoes a thermonuclear runaway when 
approaching Mch and calculate synthetic observables to compare 
the optical display of such an explosion to observed SNe Ia(x). 
For this purpose we neglect the uncertainties regarding the exact 
structure of the progenitor and follow a simplified approach to 
model the hybrid CONe WD. Motivated by the model of [Denis-] 


[senkov et al.[ 1 2013[ l, we choose the initial WD (central density 
Pc = 2.9 X 10^ g cm“^ and isothermal temperature structure with 
T = 5 X 10® K) to have a C-rich core of 0.2Mq (2f(C) = 0.5, 
X(0) = 0.5), surrounded by an ONe layer (X(C) = 0.03, 
X(0) = 0.5, X(Ne) = 0.47) up to 1.1 Mq. On top of this, we 
place a layer of accreted material with XjC) = 0.5, Tf (O) = 0.5 
which extends up to 1.4 M©. 


3 EXPLOSION SIMULATION 

To simulate the explosion phase we use our 3D hydro code LEAFS 
(based on the implementation by [Reinecke et al.|[2002[ l. LEAFS 
solves the reactive Euler equations and models the propagation of 
deflagrations using a levelset scheme i Reinecke et al.|1999[[Osh^ 


[& Sethian|1988[|Smiljanovski et aF 

bulence model fSchmidt et al.|2006; 


1997[l and sub-grid scale tur- 


y to account for flame accel¬ 


eration by buoyancy- and shear-instability induced turbulent fluid 
flows. Self-gravity is dealt with by a monopole gravity solver and 
an equation of state appropriate for WD matter is used. To follow 
the evolution of the explosion ejecta to the phase of homologous ex- 
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pansion, we use the expanding grid technique described by|Rdpk^ 
( |2005^ and |R6pke et al.| ( |2006^ . For more information on the hy¬ 
drodynamics code and details concerning the burning treatment see 
IFink et al.H20T4) . 

In our progenitor WD as described in Section]^ we ignite a 
deflagration near the centre using the same asymmetric 5-kernel ig¬ 
nition configuration as in the N5def model of| Kromer et al.|j2013 ^ 
and [Fink et al.| ( [2014^ . This model of the deflagration of a near- 
Mch CO WD produces observables in good agreement with the 
Type lax SN 2005hk ( [Kromer et al.||2013) . Here, we follow the 
flame propagation in our hybrid CONe WD and the subsequent hy- 
drodynamical evolution of the ejecta up to the phase of homologous 
expansion (see Figure[TJ. 

After the burning starts from the 5 ignition kernels, the flames 
soon merge and a one-sided deflagration plume develops within 
the C-rich core. In the ONe layer above 0.2 Mq, it is assumed that 
the C deflagration flame cannot propagate because of the small C 
fraction (X (C) = 0.03). Thus, burning is switched off in the region 
initially consisting of mainly ONe matter and further out. Burning 
ceases at ~ 1.5 s past explosion. The hot ashes rise to the surface 
and wrap around the unbumt parts of the progenitor WD. Due to 
the small nuclear energy released in the burning (Enuc = 1.0 x 
10®° erg), only 0.014 Mq of material is ejected while a massive 
bound remnant of 1.39 Mq is left behind (the mass of the bound 
remnant and ejecta are determined as described by [Kromer et al.| 
|2013[ >. We stop the simulation at 100 s after ignition. At that time 
the ejecta are in homologous expansion with an asymptotic kinetic 
energy of fCkin = 1.8 x 10"^® erg. 

In a postprocessing calculation with a 384-isotope nuclear 
network ( [Travaglio et al.|[2004[ > we determine the detailed chem¬ 
ical composition of the ejecta. To this end, about 10° Lagrangian 
tracer particles are passively advected in the hydrodynamic explo¬ 
sion simulation and thermodynamic trajectories are recorded for 
each tracer particle. The tracers are placed in the WD with variable 
masses according to the prescription by [Seitenzahl et aT] ( |2010| l. 
The initial composition of the tracer particles is drawn from the 
abundance profile of our progenitor WD. 

Within the ejecta, we find 6.41 x 10“® Mq of iron-group ele¬ 
ments (IGEs), of which 3.40 x 10“® Mq are ®°Ni. The most abun¬ 
dant other species are O (3.49 x 10“® Mq), C (2.73 x 10“® Mq), 
Si (8.02 X 10“^ Mq) and S (3.56 x IQ-^ Mq). The detailed com¬ 
position of the ejecta is given in Tables[^an d|2] 

Using an SPH-like algorithm ( [Kromer et al.|2010l ) we map the 
final composition of the unbound tracer particles to a 200® Carte¬ 
sian grid to determine the abundance structure of the ejecta (see 
Figure]^. Our model shows significant deviations from spherical 
symmetry which is a consequence of the turbulent burning and the 
one-sided propagation of the thermonuclear flame within the pro¬ 
genitor WD. 

For the bound remnant, which is not resolved in the hy¬ 
drodynamics calculations, our nucleosynthesis post-processing is 
not reliable, since the tracers do not record accurate thermody¬ 
namic trajectories. Thus, we can only give approximate yields 
from the simplified description of nuclear reactions implemented 
in our hydro scheme. From this we find 1.33 Mq unbumt material, 
3.3 X 10“® Mq intermediate-mass elements and 5.3 x 10“® Mq 
IGEs, of which 2.0 x 10“® Mq are ®°Ni. 


Table 1. Asymptotic chemical yields in the ejecta of our model. 



Mq 


Mq 

Total 

1.40 X 10-2 

Cl 

1.00 X 10-s 

He 

3.58 X 10-2 

Ar 

6.29 X 10-° 

Li 

3.90 X 10-13 

K 

6.98 X 10-2 

Be 

4.58 X 10-1® 

Ca 

4.74 X 10-® 

B 

1.83 X 10-11 

Sc 

3.24 X 10-9 

C 

2.73 X 10-® 

Ti 

1.68 X 10-° 

N 

1.79 X 10-2 

V 

9.35 X 10-2 

0 

3.49 X 10-® 

Cr 

8.71 X 10-® 

F 

2.58 X 10-9 

Mn 

1.57 X 10-1 

Ne 

3.88 X 10-1 

Fe 

5.28 X 10-® 

Na 

1.75 X 10-° 

Co 

9.05 X 10-° 

Mg 

7.60 X 10-® 

Ni 

8.76 X 10-1 

A1 

3.78 X 10-° 

Cu 

5.74 X 10-® 

Si 

8.02 X 10-1 

Zn 

6.15 X 10-® 

P 

1.95 X 10-° 

Ga 

1.20 X 10-9 

S 

3.56 X 10-1 

Ge 

1.63 X 10-9 


Table 2. Yields of radioactive isotopes in the model ejecta alt = 100 s. 


Mq Mq 


56Ni 

3.40 

X 

10- 

-3 

59 Fe 

8.15 

X 

10- 

-8 

57Ni 

1.29 

X 

10- 

-4 

52Mn 

4.95 

X 

10- 

-8 

95 Co 

1.18 

X 

10- 

-4 

41 Ca 

4.16 

X 

10- 

-8 

52 Fe 

4.36 

X 

10- 

-5 

44Ti 

3.95 

X 

10- 

-8 

55 Fe 

3.62 

X 

10- 

-5 

60 Co 

3.02 

X 

10- 

-8 

52Co 

1.82 

X 

10- 

-5 

63Ni 

1.79 

X 

10- 

-8 

53 Fe 

7.84 

X 

10- 

-6 

26 A1 

1.34 

X 

10- 

-8 

59Ni 

7.56 

X 

10- 

-6 

35 s 

2.73 

X 

10- 

-9 

53 Mn 

3.96 

X 

10- 

-6 

65 Zn 

2.42 

X 

10- 

-9 

56 Co 

2.21 

X 

10- 

-6 

49v 

2.33 

X 

10- 

-9 

iSCr 

1.17 

X 

10- 

-6 

39 Ar 

2.22 

X 

10- 

-9 

51 Mn 

6.28 

X 

10- 

-7 

32p 

2.12 

X 

10- 

-9 

32Ar 

1.58 

X 

10- 

-7 

36 Cl 

2.00 

X 

10- 

-9 

60 Fe 

1.48 

X 

10- 

-7 

33p 

1.60 

X 

10- 

-9 

62 Zn 

1.33 

X 

10- 

-7 

48 Y 

1.00 

X 

10- 

-9 

58 Co 

1.21 

X 

10- 

-7 

40k 

6.25 

X 

10“ 

■10 

5lMn 

1.12 

X 

10- 

-7 

22Na 

3.32 

X 

10“ 

■10 

19 Cr 

1.11 

X 

10- 

-7 

®2Si 

2.81 

X 

10“ 

■10 

14c 

1.07 

X 

10- 

-7 

68 Ge 

7.98 

X 

10“ 

■ 11 

51 Cr 

9.11 

X 

10- 

-8 

65Ga 

1.43 

X 

10“ 

■ 11 


4 COMPARISON WITH OBSERVATIONS 

Using our 3D Monte Carlo radiative transfer code ARTIS ([Kromer[ 
|& Sim[2009[|Sim[2007| l we calculate synthetic observables for the 
ejecta. To this end we re-map the abundance and density structure 
as obtained from the explosion simulation to a 50® Cartesian grid. 
On this grid we follow the propagation of 1.2 x 10® Monte Carlo 
quanta for 136 logarithmic time steps from 0.3 to 35 d past ex¬ 
plosion. Escaping Monte Carlo packets are binned in time and on 
a logarithmic wavelength grid spanning 10,000 bins from 600 to 
30,000 A to obtain a spectral time sequence. For our simulation we 
use the atomic data set as described by [Gall et aT] \2Q\2) . Local 
thermodynamic equilibrium is assumed for t < 0.42 d and a grey 
approximation is used in optically thick cells (cf. [Kromer & Sim[ 
|2009[ ( to speed up the calculations. 

Synthetic broad-band light curves of our model are shown 
in Figure To examine the effects of asymmetries in the ex¬ 
plosion ejecta (cf. Figure |^, we extract light curves along 100 
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5.0xl0"cm 


4.5xl0®cm 


5.0xl0“cm 


Figure 1. Hydrodynamical evolution of our model at 0.85, 2.5 and 10 s after explosion (from left to right). We show a volume rendering of the mean atomic 
number where yellow corresponds to unburnt material, and cyan and dark blue to ejecta rich in intermediate-mass or iron-group elements, respectively. 



Figure 2. Final ejecta structure at 100 s. The left-hand panel shows a 3D volume rendering of the mean atomic number as in Figurej^jthe triad corresponds to 
a length of 5 X 10^® cm). The right-hand panel shows angle-averaged density and mass fractions of selected species as a function of expansion velocity. 


different viewing angles sampling the full solid angle. The peak 
magnitudes of our model span a range from —13.2 to —14.6 and 
— 14.2 to —14.8 in B and V band, respectively. These are signif¬ 
icantly fainter than normal SNe la or proto-typical SNe lax like 
SN 2002cx or SN 2005hk. Our peak magnitudes are, however, 
in good agreement with those of the faint Type lax SN 2008ha 
(Mmax(B) = -13.74mag, Minax(’l^) = -14.21 mag; |Foley 
|et al.|2009t . We also find reasonable agreement for the 
light colours along selected lines-of-sight. 


; maximum- 


However, our model light curves evolve significantly faster 
than SN 2008ha. This becomes obvious comparing the rise times 
of the model (fmax(B) = 2.9 ... 4.6 d depending on the viewing 
angle) and SN 2008ha (~10d; [Foley et al.|2009^ , and also the de¬ 
cline of the redder bands at epochs later than ~10d past B-band 
maximum. The fast evolution of our model is a consequence of 
the extremely low ejecta mass of 0.014 M©, which is significantly 
smaller than the ejecta mass derived for SN 2008ha. Using analyt¬ 


ical light curve models [Foley et ar] ( |2009| l inferred an ejecta mass 
of Mej = 0.15 Mq based on a first set of post-maximum obser¬ 


vations. Including a pre-maximum spectrum to the analysis Foley 


et al. 


2009 


j2010^ revised this estimate to Mej = 0.30 Mq ([Valenti et al. 


obtain an ejecta mass in the range 0.1 to 0.5 Mq). 


As a consequence of the relatively strong asymmetry of the 
model ejecta, our light curves show stronger sensitivity to different 
viewing angles than those of the N5def model ( [Kromer et al.|201^ . 
At peak the U -, B- and U-band magnitudes show a scatter of 2.3, 
1.4 and 0.6 mag, respectively (Figure [^. Given the small sample 
size of faint SNe lax and the large intrinsic diversity in the few 
systems known, this viewing angle sensitivity cannot be used to 
rule out the model. 


Selected snapshots of the synthetic spectral time series result¬ 
ing from our simulations are shown in Figure [^ and compared to 
observed spectra of SN 2008ha at corresponding epochs. At —1.2 d 
(with respect to B-band maximum) our model agrees fairly well 
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Figure 3. Synthetic light curves of our model for different filter bands as indicated in the panels. We show light curves from 100 different viewing angles 
(light red) sampling the full range of solid angle and the angle average (dark red). For comparison we show observational data for the normal SNe la 2004eo, 
2005cf (grey squares and triangles, respectively), SN 2005hk (a proto-typical SN lax, dark grey circles) and SN 2008ha (black upside-down triangles). The 
latter marks the faint end of the current sample of SNe lax. 



Rest wavelength in A 

Figure 4. Angle-averaged synthetic spectra (in red) of our model at different 
epochs (times with respect to f?-band maximum are indicated by the labels 
in the panels). For comparison, observed spectra of SN 2008ha ^Foley et al.| 
|2009| |2010| [Valenti et al.|200^ are shown (the flux calibration has been 
checked against the photometry and adjusted if necessary). The observed 
spectra have been de-reddened and de-redshifted. 


with the spectral energy distribution of SN 2008ha. However, the 
model fails to reproduce some of the narrow line features present 
in SN 2008ha, e.g. the features at 6300A, 7600A, and 8400A 
that have previously been identified with Si II, O I and Ca II, respec¬ 
tively jFoley et al.|2010^ . This could indicate a lack of intermediate- 
mass elements in the outer layers of the model ejecta. 

Another explanation for the lack of narrow features in our 
model could be too large ejecta velocities. The model spectra are 
dominated by fluorescent emission in a plethora of atomic lines of 
IGEs, which are present throughout the ejecta in our model (see 
Figure]^. As long as the emission originates from the outer lay¬ 
ers (where the expansion velocity is large) the emission lines are 
blended giving rise to a smooth pseudo continuum. With time the 
ejecta expand further, the outer layers become increasingly opti¬ 
cally thin and deeper layers at lower velocities start to dominate 
the emission. Then individual lines de-blend and narrow features 
appear. A larger ejecta mass at similar energy release could reduce 
the ejecta velocities and thus lead to narrower line features. 

A week after B-band maximum (7.5 d) the agreement between 
our model and SN 2008ha is fairly good. The model shows a for¬ 
est of narrow line features similar to that observed in SN 2008ha 
though not every feature is reproduced perfectly. However, this can¬ 
not be expected from a model that is not tuned to fit the data. Even 
for normal SNe la, “first-principle” models do not perfectly repro¬ 
duce the observed spectra (e.g. jSim et al.|20l'3| >, particularly if plot¬ 
ted in absolute fluxes as in Eigure|^ In previous work the narrow 
line features were interpreted as P-Cygni absorption profiles (e.g. 
[Foley et al.||2009| >. In our model most of the features can be at¬ 
tributed to fluorescent emission in lines of IGEs (see Figure]^. 

At later epochs (15.1 d and 22.5 d) the agreement deteriorates. 
In particular, the rapid decrease of flux at wavelengths larger than 
5500 A is not observed in SN 2008ha. This reflects the rapid decline 
in some of our broad-band light curves {R band and redder) dis¬ 
cussed above and indicates that the model ejecta become optically 
thin too early, again pointing at a lack of ejecta mass in our model 
compared to SN 2008ha. At late epochs, however, some approxi¬ 
mations on the plasma state, made in our radiative transfer code. 
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Figure 5. Angle-averaged synthetic spectrum of our model (grey) at 7.5 d 
past f?-band maximum. The colour coding shows the chemical species re¬ 
sponsible for both bound-bound emission and absorption of quanta in our 
Monte Carlo simulation. Below the synthetic spectmm, we colour code the 
fraction of escaping quanta in each wavelength bin that were last emitted by 
a particular chemical species (corresponding atomic numbers are illustrated 
in the colour bar). The coloured regions along the top of the plot indicate 
which elements were last responsible for removing quanta from a particular 
wavelength bin (either by absorption or scattering/fluorescence). 

become less applicable ( [Kromer & Sim|2009^ . Indeed, we find that 
the agreement between synthetic observables and data deteriorates 
also for other model classes (e.g. jRdpke et al.|2012^ . Thus, some of 
the differences at later epochs could also be due to systematics in 
the radiative transfer. 


5 DISCUSSION 

5.1 Ejecta mass 

The lack of narrow line features around maximum light and the 
overly rapid light curve evolution of our model compared to 
SN 2008ha indicate that the ejecta mass of the model is too low. 
In contrast, the good agreement in peak luminosity suggests that 
the ®®Ni yield in the model ejecta is close to that in SN 2008ha. 
Here, we have investigated a single realisation of a deflagration in 
a hybrid CONe WD for a particular progenitor structure and igni¬ 
tion setup. In reality both of these will show some diversity. The 
progenitor structure will depend on the accretion history and the 
final reactive-convective evolution of the WD core before the ther¬ 
monuclear runaway, which will also determine the ignition setup. 
A thorough investigation of the parameter space is necessary to as¬ 
sess whether variations in these initial conditions can increase the 
ejecta mass while keeping the ®®Ni mass in the ejecta similar, and 
thus lead to a better agreement with SN 2008ha. This requires a 
large number of explosion simulations and will be addressed in fu¬ 
ture work. Given the large uncertainty in the chemical and thermal 
structure of hybrid WDs ( jPenissenkov et al.|2015| ), here we focus 
on a simple model for an exploratory study. 

5.2 Bound remnant 

A key prediction of our model is a bound stellar remnant of ~ 
Mch- Recently, deep post-explosion HST images have revealed a 


Figure 6. Bolometric light curve of the ejecta (red) and instantaneous en¬ 
ergy deposition due to radioactive decay of ®®Ni in the bound remnant 
(blue). For comparison the green line shows the Eddington luminosity of 
a Mqu object. To obtain the bolometric light curve of the ejecta for the full 
time range, we have performed an additional low-resolution radiative trans¬ 
fer simulation from 2 to 120 d to extend our high-resolution simulation that 
stops at 35 d. 


point source S1 at the site of the faint Type lax SN 2008ha i Foley 
et al.|2014^ . One possible interpretation of SI, proposed by Foley 
et al.j ( |2014bt , is that it could be an inflated stellar remnant of a 
failed deflagration. From our present simulations, which follow the 
evolution of the explosion ejecta on an expanding grid, we cannot 
resolve the structure and thermal state of the remnant. Thus we can¬ 
not compare our bound remnant to the observed properties of the 
point source. We note, however, that explosion energy deposited in 
the bound remnant could lead to an expanded envelope. Whether 
the expansion is sufficient to account for a radius of 1500 Rq as 
derived for SI ( [Foley et al.|2014b| ( remains to be seen. This ques¬ 
tion will be addressed in more detail with an adaptive grid code in 
future work. 

Moreover, our simulations show that 2.0 x 10“^ Mq of ®®Ni 
synthesized during the explosion stay in the bound remnant (see 
Section]^. This is ~ 6 times larger than the ®®Ni yield in the ejecta. 
Consequently, the instantaneous energy deposition by decaying 
®®Ni in the bound remnant is significantly larger than the bolo¬ 
metric luminosity of the ejecta derived from our radiative transfer 
simulations (see Figure]^. Again, our ignorance of the structure of 
the bound remnant makes it difficult to make a robust statement of 
how the decay energy will affect the observational display. If the 
®®Ni is located close to the surface of the remnant, it might drive 
further mass ejection and contribute significantly to the observed 
luminosity. In contrast, if most of the ®®Ni is confined to the cen¬ 
tral regions of the stellar remnant, it might take a long time for the 
decay-deposited energy to diffuse to the surface of the remnant so 
that it becomes negligible for the observational display of the ex¬ 
plosion. 

Another interesting question concerns the long-term evolution 
of the bound remnant. Because the kinetic energy of the explosion 
is very low, the donor star is presumably not much affected by the 
impact of the ejecta ( [Liu et al.|20T3) . Thus, it is possible that ac¬ 
cretion resumes once the remnant has relaxed into an equilibrium 
state. Since only a tiny fraction of the mass of the initial WD is 
ejected, it is likely that the bound remnant approaches the Mch 
limit again. Whether or not this may lead to subsequent deflagra¬ 
tions and thus a recurrent nature of some SNe lax, depends on the 
chemical and thermal structure of the bound remnant. If the C frac- 
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tion in the bound remnant is low, an accretion-induced collapse and 
a neutron star might be the more probable outcome. As our current 
simulations are unable to resolve the bound remnant, we can not 
make a predictive statement about its final fate at this stage. 


5.3 Burning in the ONe layer 

In this work, we have assumed that deflagration burning ceases 
when the flame front reaches the ONe mantle of the hybrid WD. 
This assumption significantly reduces the energy release compared 
to deflagrations in CO WDs and is critical to obtain very low ®®Ni 
masses as observed in faint SNe lax. However, in principle, defla¬ 
grations are possible in ONe material as well jTimme s & Woosley| 
[T992l >. Since the energy release from burning this fuel is lower, 
the flame slows down and its width increases more rapidly with 
lower fuel density than for deflagrations in CO material jTimmes| 
|& Woosley|1992[ l. ONe deflagrations are therefore not expected to 
propagate at densities lower than ~ 10® g cm“^. In our model, the 
CO deflagration reaches the ONe layer when it has expanded to 
densities of ~ 1.2 x 10® g cm“®, so some additional burning in the 
ONe material is possible. Detailed microscopic flame simulations 
are necessary, however, to assess our assumption that burning stalls 
shortly after the deflagration reaches the ONe layer. 


5.4 Binary population synthesis - rates and constraints on 
delay times 

It is important to assess the likelihood of such events from a the¬ 
oretical standpoint: how frequent are the faint SN lax events? We 
evolved 300,000 binaries with Z=0.02 from the ZAMS up to a Hub¬ 
ble time assuming a binary fraction of 70 percent using the popula¬ 
tion synthesis code StarTrack ( [Belczynski et al.|2002||200^ . To 
obtain theoretical birthrates we first calculated the number of ONe 
WDs that approach Mch via stable Roche-lobe overflow from a 
stellar companion (see P-MDS model description in |Ruiter et al.| 
|2014| l. Typically, these systems are considered to lead to accretion- 
induced collapse and form neutron stars, but as discussed in Sec¬ 
tion]^ if these WDs contain some non-negligible fraction of un- 
bumt C, they may instead lead to thermonuclear explosions. 

The StarTrack code currently does not account for the evo¬ 
lution of hybrid CONe WDs explicitly. However, it is reasonable 
to assume that if such hybrid WDs exist then the lower mass limit 
will occur near the boundary where, in our population models, a 
degenerate CO core is formed, and where C burning occurs non- 
explosively leading to the formation of a degenerate ONe core (the 
CO WD - ONe WD boundary; see [Belczynski et ar]j2008| l). The 
upper limit for the hybrid core mass will lie somewhere within the 
range of masses that are canonically assumed to result in ‘pure’ 
ONe WDs. 

jPenissenkov et al^ ( |2015^ found that single stars with ZAMS 
masses between ~ 6.4 and 7.3 M© produce CONe hybrid WDs. 
This same mass range cannot be extrapolated to interacting binary 
stars since a star that has lost or gained mass will follow a different 
course of evolution (and end up with a different core mass) than that 
of a single star with the same ZAMS mass. To estimate how many 
of our ONe WDs may contain some fraction of unbumt C in their 
cores, we checked the corresponding WD birth masses that arise 
from ZAMS single stars within this mass range in StarTrack. 
The corresponding range is 1.193 to 1.325 Mq. Here, we assume 
these (ONe-rich) WDs contain some fraction of unburnt C and thus 
are hybrid WDs. A small number of ONe WDs are found below this 



Figure 7. Delay time distribution of CONe WDs that approach Mqu 
due to accretion from a binary star companion. Blue systems are those 
with helium-burning stars as donors while the red systems contain main- 
sequence, sub-giant, giant or WD donors. Numbers (y-axis) are not scaled 
to an absolute rate but rather represent the original numbers from our model. 
An estimate of absolute rates (over a Hubble time) is given in the text. 


mass boundary in our model, and so we include these as potential 
hybrid cores as well. We assume all of these WDs undergo an off- 
centre deflagration once they approach Men- 

In terms of relative rates for different SN la progenitors, we 
find that over a Hubble time, hybrid CONe WDs that may produce 
faint lax-like events are 1 percent of the rate of the entire CO-CO 
WD merger population. By comparison, they have about the same 
relative rate that we predict for the classic single-degenerate sce¬ 
nario at near solar metallicity, whereby a CO WD accretes toward 
Mch from a hydrogen-burning star (cf. table 1 of Marquardt et ah, 
submit ted). To put it in a more absolute context: [Badenes & Maozj 
1 2012 I quote a SN la rate of 1.1 x 10“^® yr“^ M©”^ for Milky 
Way like galaxies. We find from our population synthesis model 
that the CO-CO WD merger rate (averaged over a Hubble time) is 
1.06x 10“^® yr“^ Mq“^, mother words: very close to the Sbe-like 
galaxy SN la rate (see also [Li et al.|2011a^ . Taking this number at 
face value as the overall SN la rate, we find the rate of deflagrations 
in hybrid CONe WDs to be on the order of 1 percent of the SN la 
rate. This relative rate will increase, if one considers galaxies with 
active star formation rather than older stellar populations (like the 
Milky Way). Given the large uncertainties in the observed rate of 
SNe lax (different authors give values between 5 and 30 percent of 
the overall rate of SNe Ia,|Li et al.|201 lb||Foley et al.|2013||White| 
jet al.|2015f, our estimated rate seems in rough agreement with faint 
SNe lax. 

In Figure]^ we show the delay time distribution for the pop¬ 
ulation of hybrid CONe WDs estimated from our population syn¬ 
thesis model. One third of our hybrid systems have delay times 
< 150 Myr with the shortest delay time occurring at 30Myr. All 
of our prompt systems below 150 Myr have helium-burning stars 
as donors (shown in blue). The range of delay times for these pro¬ 
genitors agrees with the results of [Wang et al.jpOlT^ , who esti¬ 
mated the delay time range for CONe hybrid WDs that accrete to¬ 
ward Mch from helium stars to be 28 to 178 Myr. Our results are 
consistent with the fact that SNe lax are found among young stel¬ 
lar populations; SN 2008ha is estimated to have a delay time of 
80Myr (Foley et al. 2014b i. Nearly half of our hybrid systems, 
however, have hydrogen-burning donors, while 17 percent have He 
WD donors (all shown in red). All of these systems have longer 
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delay times (shortest one is 420 Myr). The main factor that deter¬ 
mines whether a progenitor will have a very short delay time and a 
helium-huming star donor, or whether it will have a slightly longer 
delay time with a hydrogen-burning star (or WD) donor is the ini¬ 
tial mass of the donor star. For the prompt systems, the ZAMS mass 
of the donor star ranges from 5 — 10 Mq while for the rest, it ranges 
from 1 — 3 Mq. Given our findings, it is unlikely that most progen¬ 
itors of faint SN lax events originate from systems with hydrogen- 
burning or WD donors, though our results strongly support the idea 
that faint SN lax events arise from CONe WDs that accrete from 
helium burning stars. 


6 CONCLUSION 


Previous models of deflagrations in near-Mch CO WDs had dif¬ 
ficulties in explaining faint SNe lax, since they produced signifi¬ 


cantly too large amounts of ®®N1 iFink et al.|2014 


_ _ . To reduce the 

®®Ni production, burning has to quench at an early stage. Defla¬ 
grations in hybrid CONe WDs, which have recently been predicted 
by stellar evolution models (|Denissenkov et al.|2013| |Chen et al.| 
[ 20 T^ , may lead to such a quenching when the burning front tran¬ 
sitions from the WD’s CO-rich core to its ONe-rich mantle. 

We have simulated an off-centre deflagration in a near-Mch 
hybrid CONe WD. Assuming that the deflagration cannot propa¬ 
gate in the ONe layer, we find fhat only 0.014 Mq of material is 
ejected while the remainder of the mass stays bound. The ejecta 
consist predominantly of IGEs, O, C, Si and S. Containing only 
3.4 X 10“® Mq of ®®Ni, the ejecta will give rise to a faint tran¬ 
sient. 

Performing radiative transfer simulations for our model ejecta, 
we find peak absolute magnitudes in the range M^ax = —13.2 
to —14.6 mag depending on the viewing angle. This is in good 
agreement with the observed peak brightness of the faint Type lax 
SN 2008ha (M®ax = —13.7mag). Between peak and 10 d there¬ 
after we also find reasonable agreement between the observed spec¬ 
tral shape of SN 2008ha and our model. This shows for the first time 
that deflagrations in near-Mch hybrid CONe WDs can lead to faint 
transients with a display similar to faint SNe lax. 

However, our model still has some shortcomings. At epochs 
before maximum and later than 10 d past maximum the spectral 
evolution is too fast, probably indicating that the ejecta mass is too 
low. Since the exact thermal and chemical structure of Mch hybrid 
CONe WDs is not well constrained jPenissenkov et al.|2015t , we 
used a simple model for the progenitor. Future work exploring a 
range of different initial conditions will show whether even better 
agreement with faint SNe lax is possible from deflagrations in hy¬ 
brid CONe WDs. More vigorous ignition setups, similar to those 
explored in CO WDs by [Fink et aL] ^2014^ , or hybrid WDs with 
more massive CO cores might also be able to explain intermediate- 
luminosity SNe lax like SN 2007qd ( [McClelland et al.|2010| ). 

From binary population synthesis calculations we find that 
near-Mch hybrid CONe WDs that may produce faint SNe lax oc¬ 
cur at a rate of 1 percent of the Galactic SN la rate. The delay 
times of these systems depend sensitively on the nature of the donor 
star. He-buming donors (delay times 30 to 150 Myr) are compatible 
with the estimated delay time of SN 2008ha (80 Myr, [Foley et ah] 
|2014a [ and the observed preference for star-forming environment 
among SNe lax ( [Lyman et al.||201^ . H-burning and WD donors 
have somewhat longer delay times. 

Finally, hydrodynamic simulations with adaptive grid codes 
will be needed to resolve the structure and thermal state of the 


bound remnant predicted by our model. This will be crucial to as¬ 
sess the influence of fhe bound remnant on the observational dis¬ 
play of the explosion and allow for a detailed comparison to the 
potential stellar remnant detected at the site of the faint SN 2008ha. 
Following the long-term evolution of the near-Mch bound remnant 
and companion star will show whether subsequent mass-transfer 
episodes can lead to recurrent explosions or accretion-induced col¬ 
lapse. 
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